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Abstract. We study the predictions for structure formation in an induced gravity dark energy 
model with a quartic potential. By developing a dedicated Einstein-Boltzmann code, we study 
self-consistently the dynamics of homogeneous cosmology and of linear perturbations without 
using any parametrization. By evolving linear perturbations with initial conditions in the 
radiation era, we accurately recover the quasi-static analytic approximation in the matter 
dominated era. We use PLANCK 2013 data and a compilation of baryonic acoustic oscillation 
(BAO) data to constrain the coupling 7 to the Ricci curvature and the other cosmological 
parameters. By connecting the gravitational constant in the Einstein equation to the one 
measured in a Cavendish-like experiment, we find 7 < 0.0012 at 95% CL with PLANCK 2013 
and BAO data. This is the tightest cosmological constraint on 7 and on the corresponding 
derived post-Newtonian parameters. Because of a degeneracy between 7 and the Hubble 
constant Hq, we show how larger values for 7 are allowed, but not preferred at a significant 
statistical level, when local measurements of Hq are combined in the analysis with PLANCK 
2013 data. 
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1 Introduction 


Inflation or quintessence are naturally embedded in scalar-tensor theories of gravity. In these 
models the scalar field which regulates the gravitational coupling also drives the acceleration of 
the Universe. The non-minimal coupling to gravity can change significantly the perspective 
on inflation or quintessence in Einstein general relativity. In the inflationary context, for 
instance, a large coupling of the inflaton to gravity allows potentials with a self-coupling 
which would be excluded in the minimally coupled case [1-4]. In the dark energy context, for 
example the coupling to gravity could allow super-acceleration (i.e., H > 0) with standard 
kinetic terms for the scalar held [5]. 

In this paper we consider induced gravity (IG) with a quartic potential V(a) = Acj 4 /4 
as a simple scalar-tensor dark energy model: 
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where C m denotes the contribution by matter and radiation. Under a simple held redefinition 
7 <r 2 = 0 /( 87 t), the above action can be cast in a Brans-Dicke-like model [ 6 ] with a quadratic 
potential: 
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with the following relation between the dimensionless parameters of the two theories: 


1 \/2A 

“ BD = ^ 


(1.3) 


The action in Eq. (1.1), which contains only dimensionless parameters, was introduced to 
generate the gravitational constant and inflation by spontaneous breaking of scale invariance 
in absence of matter [7]. In the context of late cosmology, this action was studied in Refs. [ 8 , 9] 
to reduce the time dependence of the effective gravitational constant in the original Brans- 
Dicke model (i.e., with a vanishing potential [ 6 ]) and to generate an effective cosmological 
constant. The cosmological background dynamics from Eq. (1.1) was shown to be consistent 
with observations for small 7 , i.e., 7 < 1CU 2 [10]. 
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The potential term in Eq. (1.1) is important for the global dynamics of the model 
and modifies the original Brans-Dicke attractor with power-law time dependence of the 
scalar held in presence of non-relativistic matter, i.e., a(t) = (t/to)( 2 ‘ JBD+2 )/( 3a -’ BD + 4 ) and 
= < M^Ao) 2/(3 ‘ Jbd+4) - At recent times, the potential term drives the Universe into acceler¬ 
ation and Einstein gravity plus a cosmological constant with a time-independent value of the 
scalar held emerge as an attractor among homogeneous cosmologies for the model in Eq. (1.1). 

In this paper we study structure formation in the IG dark energy with a quartic potential 
in Eq. (1.1). We study how gravitational instability at linear level depends on 7 through a 
dedicated Einstein-Boltzmann code. We then use these theoretical predictions for cosmolog¬ 
ical observables to constrain the model with the PLANCK 2013 data [11—13], a compilation 
of baryonic acoustic oscillations (BAO) data [14-16] and local measurements of the Hubble 
constant [17-19]. 


2 Dark Energy within Induced Gravity 


The Friedmann and the Klein-Gordon equations for IG in a hat Robertson-Walker metric are 
respectively: 


a 3y < 7 2 670- 2 
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once the Einstein trace equation: 
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is used. In the above V a denotes the derivative of the potential V(cr) with respect to <7, 
the index i runs over all huid components, i.e. baryons, cold dark matter (CDM), photons 
and neutrinos, and we use a dot for the derivative with respect to the cosmic time. When 
considering V oc c 4 the potential cancels out from the Klein-Gordon equation and the scalar 
held is driven by non-relativistic matter. In the rest of the paper we will restrict ourselves to 
V(a) = Aa 4 /4. 

We consider the scalar held a at rest deep in the radiation era, since an initial non¬ 
vanishing time derivative would be rapidly dissipated [10]. The scalar held is then driven by 
non-relativistic matter to an asymptotically value higher than the one it had in the radiation 
era as can be seen in the left panel of Fig. (1); when the scalar held freezes the Universe is 
driven in a de Sitter era by the scalar held potential which behaves as an effective cosmo¬ 
logical constant [ 8 , 10], as can be seen in the central panel of Fig. (1). Since a regulates 
the gravitational strength in the Friedmann equations, the present value of the held 07 can 
be connected the gravitational constant G measured in laboratory Cavendish-type and solar 
system experiments by the relation: 


2 _ 1 I + 87 

7f7 ° 87 tG 1 + 67 


(2.4) 


where G = 6.67 x 10 ~ 8 N cm 3 g _1 s~ 2 . The above equation assumes that er is effectively 
massless on Solar System scales. Note that the scalar held a is effectively massless in the 
radiation dominated era, as can be seen in Eq. (2.2). 
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Figure 1. Evolution of er/cro (left panel), fij (middle panel) and wde (right panel) as function of 
ln(a) for different choices of 7 . 


The evolution of the background cosmology can be easily compared with dark energy 
in Einstein gravity with a Newton’s constant G]\r given by the scalar field value at present 
G = (87T7<7o) _1 . The Friedmann equation can be therefore rewritten by introducing an 
effective dark energy component [20], whose energy and pressure densities for this model are 
[ 10 ]: 
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In the central panel of Fig. (1) we display the time evolution of the density contrasts of 
radiation - Or = (p u + p 7 )/(37<7Q.£f 2 ) -, matter - 17m = (Pb + Pcdm) /(S'l&oH 2 ) - and effective 
dark energy - Ode = PDE/i^'l^oH 2 )) As a third panel in Fig. (1), we display the time 
evolution of the parameter of state of the effective dark energy component, u>de = Pde/pde- 


3 The evolution of cosmological fluctuations 


As for the background, we study linear fluctuations in the Jordan frame. We consider metric 
fluctuation in the longitudinal gauge: 


ds 2 = —dt 2 ( 1 + 2 T(t, x)) + a 2 (t){ 1 — 2 4>(f, ~x))dx l dxi (3.1) 


and for the scalar held: 


cr(t , x) = a(t) + 5a(t, x). 
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The perturbed Einstein equations for our IG model with a quartic potential in the longitudinal 
gauge are: 
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In the above pi,pi ( 5pi,5pi ) denote the energy and (longitudinal) pressure density pertur¬ 
bations for each matter component, respectively. The velocity potential and the anisotropic 
stress are denoted by Oi and cf,;. We refer to Ref. [21] for the conservation of the CDM, baryons, 
photons and neutrino energy-momentum tensors, since these equations are unchanged from 
those in Einstein gravity. 

The Klein-Gordon equation at linear order in the longitudinal gauge is: 
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It is interesting to note that the equation for the field fluctuation does not depend on 
the potential explicitly for the self-interacting case, as for the background in Eq. (2.2). 

We have modified the publicly available Einstein-Boltzmann code CLASS 1 [22, 23] to 
evolve background and linear fluctuations within induced gravity. Previous implementations 
of induced gravity in Einstein-Boltzmann codes include Refs. [24-28]. 

We initialize the fluctuation of the metric and of the matter components with adiabatic 
initial condition deep in the radiation era. We have tested our numerical results from our 
modified code against analytic approximations derived within the matter era. To this purpose 
we consider the quantity p(k, a) which parametrize the deviations of \L from Einstein gravity 
[29, 30]. We consider the definition for p(k, a) which holds also during the radiation dominated 
regime as in Ref. [31]: 

A 2 ’E = —AirGa 2 p(k, a) [A + 3 (p + p)a\ (3.5) 

where A = Yli^Pi + 3 aH(pi + pi)0i/k 2 , with 5i = 5pi/pi and 6i is the velocity potential. 
Analogously we consider the deviations from Einstein gravity of the difference between the 
Newtonian potentials, parametrized by 5, whose definition valid also in the radiation domi¬ 
nated regime is [31]: 

k 2 [& — 5(k, a)Af] = 12irGa 2 p(k, a)(p + p)a (3.6) 

Our results are shown in Fig. (2). In the left panel we show the evolution of p(k,a) 
for two wavenumbers (k = 0.05,0.005 Mpc -1 ) and two values of the coupling to the Ricci 

1 www. class-code.net 
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Figure 2. Comparison of theoretical approximations for /i and <5 which parametrize deviations from 
Einstein gravity (black lines) with our numerical results for two wavenumbers (fc Mpc = 0.05,0.005) 
and two values of the coupling to the Ricci curvature (7 = 10 ~ 2 ,10 -3 ). 


curvature (7 = 10 2 ,10 3 ). We compare our numerical results for fi(k,a) to the analytic 
approximation in the matter era: 

H(k,a) = - 1 (3.7) 

which is derived from Ref. [30] for our choice of the potential and for our identification of the 
gravitational constant in Eq. (2.4). Well after matter-radiation equivalence, the quasi-static 
analytic approximation for k 3> aH in the matter era for fi(k. a) independent on k is well 
recovered. In the right panel of Fig. (2), we compare the evolution of S(k, a) for the same two 
wavenumbers and two values of the coupling 7 with the quasi-static approximation [5, 32, 33]: 


6(k, a) 


I + 47 
I + 87 


(3.8) 


Again, the analytic quasi-static approximation holds well after matter-radiation equality for 
sub-Hubble scales. 

The agreement between our numerical treatment and the quasi-static approximation 
means that our self-consistent treatment of background and linear perturbations is sufficiently 
ready for precision cosmology. The two panels in Fig. (2) also show how the time evolution for 
H and 5, independent on k, recovered within the quasi-static approximation is not valid when 
the wavelength is larger than the Hubble radius: predictions for CMB anisotropies in this 
model would be therefore affected by considering 5 constant and equal to the value obtained 
within the quasi-static approximation at all times. 


4 CMB anisotropies and Matter Power spectrum 

In the left panel of Fig. (3) are shown the power spectra of the CMB temperature anisotropies 
for different values of 7 (10 -2 , 10 -3 , 10 -4 ). The relative differences with respect to the ACDM 
reference model are shown in the right panel of Fig. (3). The change in the matter-radiation 
equality present in this scalar-tensor model [34] induces relative differences in the temperature 
power spectrum at few percent level for 7 = 10 ~ 3 . 
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Figure 3. CMB temperature anisotropies power spectrum for different values of 7 (left panel) and 
relative differences with respect to a reference ACDM (right panel). 
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Figure 4. Lensing power spectrum for different values of 7 (left panel) and relative differences with 
respect to a reference ACDM (right panel). 


In Figs. (4, 5) we display the predictions for the spectrum of lensing potential and 
its correlation with the temperature held. In Fig. ( 6 ) we display the (linear) matter power 
spectrum at z = 0 and the relative differences with respect to the ACDM reference model. 
Overall, differences at the percent level are obtained for 7 = 10 -3 in different cosmological 
observables. 

5 Constraints from cosmological observations 

We explore the parameter space by the Monte Carlo code for Cosmological Parameter ex¬ 
traction Monte Python 2 [35] connected to the modified version of the Einstein-Boltzmann 
code CLASS used in the previous sections. We use the nominal mission data release from 
PLANCK, available from the Planck Legacy Archive 3 [11], The PLANCK likelihood covering 

2 www.montepython.net 

3 pIa.esac.esa.int/pIa/aio/pIanckProducts.html 
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Figure 5. Temperature-lensing cross-correlation power spectrum for different values of 7 (left panel) 
and differences with respect to a reference ACDM (right panel). 
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Figure 6. Linear matter power spectrum (at z = 0) for different values of 7 (left panel) and relative 
differences with respect to a reference ACDM (right panel). 


temperature anisotropies from t = 2 to 2500 is combined with the low-f? WMAP polarization 
data [36-38] (this combination is denoted as Planck 2013 in the following). 

We use the PLANCK 2013 likelihood in combination with constraints on Dy(z)/r s (the 
ratio between the spherically averaged distance scale Dy to the effective survey redshift, z. 
and the sound horizon, r s ) inferred from a compilation of BAO data. These are 6dFGRS data 
[14] at z = 0.106, the SDSS-MGS data [15] at z = 0.15, and the SDSS-DR11 CMASS and 
LOWZ data [16] at redshifts z = 0.57 and 0.32. 

We vary the parameters of the flat ACDM model, i.e. the baryon density (Db/ 1 2 ), the 
CDM density (Dcdm)j the reduced Hubble parameter (h = Ho/(100 km s _1 Mpc^ 1 )), the 
reionization optical depth r, the amplitude and tilt of the primordial spectrum of curvature 
perturbations (M s and n s ) at the pivot scale k* = 0.05 Mpc" 1 . The IG dark energy model 
with quartic potential is described by these six plus one extra parameter which quantifies the 
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Figure 7. Comparison of marginalized joint 68% and 95% CL for ( r ypN,Ho) (left panel) and 
( 7 PAr,fl m ) (right panel) for Planck 2013 + Hq (grey contours), Planck 2013 + (red contours) 
and Planck 2013 + BAO (blue contours). 


coupling to the Ricci curvature 4 5 . Following Ref. [39] we sample on the quantity (, defined as: 

( = ln(l + 4j) = In ( H-— ) (5.1) 

V ^bd / 

with the prior [0,0.039] used in Ref. [39]. In this paper we consider three massless neutrinos 
’. Nuisance parameters for foreground, calibration and beam uncertainties [11, 13]. 

Our results with PLANCK 2013 + BAO data for the main and derived parameters are 
summarized and compared with the ACDM values in Table 1. The induced gravity model 
with a quartic potential is not preferred over Einstein gravity with A (Ay 2 ~ —21n£ = 0.7). 

We quote the following PLANCK 2013 + BAO 95% CL constraint on the coupling to the 
Ricci curvature: 

7 < 0.0012 (95 % CL, Planck 2013 + BAO). (5.2) 

We quote as a derived parameter the corresponding constraint on the post-Newtonian pa¬ 
rameter 7 pn = (1 + 47 )/(I + 87 ) 6 : 

0.9953 < 7pn < 1 (95 % CL, Planck 2013 + BAO). (5.3) 

It is also useful to quote the derived constraints on the change of the Newton constant 
between the radiation era and the present time 5Gn/Gn = (a 2 — (Jq)/^: 

= —0.015^q QQg (95% CL, Planck 2013 + BAO) (5.4) 

Gn 

4 The parameter of the Lagrangian A and the initial value of the scalar field at deep inside the radiation 
era are chosen to reproduce the present value of h and of the field in Eq. (2.4) by evolving the Friedmann and 
Klein-Gordon background equations. 

5 Note that the Planck collaboration assumes one massive neutrinos with a mass of 0.06 eV [13]. Given 
the interest in neutrino masses within modified gravity (see for example [40]), we will study this issue in the 
context of induced gravity in a separate publication. Even if the assumption of a mass of 0.06 eV has a small 
effect on the cosmological parameters at the Planck precision (as a 0.5 a shift to smaller value for H 0 [13]), 
we quote the results for a ACDM cosmology with three massless neutrinos in the following for a consistent 
comparison with the class of dark energy models studied here. 

6 In this class of models /?pn = 1. 
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Table 1. Constraints on main and derived parameters (at 68 % CL if not otherwise stated), 
and the constraint on its derivative ( Gn/Gn = — 2gq/(Jq) at present time: 

< ^ L {z = 0) = —0.61^0 25 [xlO -13 yr ~ l ], (95% CL, Planck 2013 + BAO). (5.5) 

Gn 

The constraints derived here are tighter than those obtained in the literature with 
PLANCK 2013 data for similar scalar-tensor models with a power-law potential [39, 41] (see 
Refs. [42-44] for analysis with pre-PLANCK data). Avilez and Skordis [41] considered the case 
of a constant potential in Brans-Dicke-like theory and quote (1 + 67 )/(l + 87 ) = 1.071q'}q at 
95% CL as the tightest constraint with a prior wbd > —3/2; we obtain [0.998,1] as the 95% CL 
range for the same quantity with PLANCK 2013 + BAO by varying £ in the interval [0,0.039]. 
Li et al. [39] considered the case of a linear potential in Brans-Dicke (i.e., a quadratic poten¬ 
tial in induced gravity) and quote 0 < £ < 0.549 x 10 -2 at 95% CL and Gat /Gjy = —1.42+2 27 
at 68 % CL from PLANCK 2013 with the same prior on £, although in combination with a 
different compilation of BAO data [39]. Note that for power-law potentials different from the 
quartic case studied here, Einstein gravity plus a cosmological constant with a independent 
on time is not the attractor at future times [45]. We therefore expect that the models studied 
in Refs. [39, 41] differ from the case of a quartic potential, in particular at recent redshifts. 

5.1 Combination with local measurements 

As from Table 1, the model considered here prefers a higher value of the Hubble parameter 
Hq with respect to ACDM. We therefore analyze the combination of the local measurements 
of the Hubble constant with PLANCK 2013 and BAO data. The local estimates of Hq are 
consistently higher than those from CMB (and BAO) and this discrepancy became a 2.5 
a tension after the PLANCK 2013 release [13]. This tension might be sign of new physics, 
although reanalysis subsequent to the PLANCK 2013 release have highlighted how hidden 
systematics and underestimated uncertainties could hide in the local measurements of Hq 
[18, 19]. For these reasons we consider separately the impact of two different local estimates 
of Hq , such as Hq = 73.8 ± 2.4 km s -1 Mpc -1 [17], denoted as 77g, and Hq = 70.6 ± 3.0 km 
s -1 Mpc -1 [19], denoted as H q. Our results are summarized in Table 2 and Fig. (7). 
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Table 2. Constraints on main and derived parameters at 68 % CL (if not otherwise stated). 


With the higher local estimate of Hq [17] we obtain a posterior on £ which is different 
at 2 a level from Einstein gravity. With the lower estimate for Hq obtained by Efstahiou 
[19] with the new revised geometric maser distance to NGC 4258 [18], we obtain a posterior 
probability for £ compatible with Einstein gravity. As can be seen by comparing the last 
columns of Table 1 and 2, the lower local estimate of Hq has almost a negligible impact when 
PLANCK 2013 and BAO data are combined. 


6 Conclusions 

We have studied structure formation in a induced gravity dark energy model with a quartic 
potential. In this model the current acceleration stage of the Universe and an accompanying 
change in the gravitational constant on large scales are due to a change in the background 
scalar held triggered by the onset of the matter dominated stage [ 8 , 10 ]. 

We have shown that the model approaches Einstein gravity plus a cosmological constant 
in the limit 7 —> 0 also at linear level. We have shown how the quasi-static parametrization 
with fj, and 7 independent on k holds only well after matter-radiation equality for sub-Hubble 
scales. 

We have derived CMB and BAO combined constraints on 7 for the case of induced 
gravity with a quartic potential. By using PLANCK 2013 [12] and BAO [14-16] data we derive 
the 95%CL constraint 7 < 0.0012, which is tighter than previous cosmological constraints on 
similar models [39, 41]. This cosmological constraint is compatible, but weaker than those 
within the Solar System [10] which can be derived by Cassini data [46]. Since there is a 
positive correlation between 7 and Ho, the combination of local measurements of Ho [17, 19] 
allows larger values of 7 , but not at statistical significant level. This analysis shows how 
a self-consistent variation of G from the radiation era to the present time can be tightly 
constrained from PLANCK 2013 and BAO data at the percent level. It will be interesting to 
see how Planck 2015 data [47] change these constraints [48]. 
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